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ABSTRACT 

The optical light-curves of GRB afterglows 990123 and 021211 exhibit a steep 
decay at 100-600 seconds after the burst, the decay becoming slower after about 10 
minutes. We investigate two scenarios for the fast decaying early optical emission of 
these GRB afterglows. In the leveise-foiwaTd shock scenario, this emission arises in 
the reverse shock crossing the GRB ejecta, the mitigation of the light-curve decay 
occurring when the forward shock emission overtakes that from the reverse shock. 
Both a homogeneous and wind-like circumburst medium are considered. In the wind- 
bubble scenario, the steeply decaying, early optical emission arises from the forward 
shock interacting with a r~^ bubble, with a negligible contribution from the reverse 
shock, the slower decay starting when the blast wave reaches the bubble termination 
shock and enters a homogeneous region of the circumburst medium. 

We determine the shock microphysical parameters, ejecta kinetic energy, and cir- 
cumburst density which accommodate the radio and optical measurements of the GRB 
afterglows 990123 and 021211. We find that, for a homogeneous medium, the radio 
and optical emissions of the afterglow 990123 can be accommodated by the reverse- 
forward shock scenario if the microphysical parameters behind the two shocks differ 
substantially. A wind-like circumburst medium also allows the reverse-forward shocks 
scenario to account for the radio and optical properties of the afterglows 990123 and 
021211, but the required wind densities are at least 10 times smaller than those of 
Galactic Wolf-Rayet stars. The wind-bubble scenario requires a variation of the micro- 
physical parameters when the afterglow fireball reaches the wind termination shock, 
which seems a contrived feature. 



1 INTRODUCTION 



There are currently two GRB afterglows for which a fast 
falling-ofF optical emission was detected at early times, only 
~ 100 seconds after the burst. The general consensus is 
that this emission arises from the GRB ejecta which is en- 
ergized by the reverse shock (RS) crossing the ejecta and 
caused by the interaction of the ejecta with the circum- 
burst medium (CBM). This interaction also drives a for- 
ward shock (FS) energizing the swept-up CBM, to which the 
later afterglow emission is attributed (the "reverse-forward 
shock" scenario). 

The RS emission was first calculated by Meszaros & 
Rees (1997), who considered the cases of a frozen- in and 
turbulent magnetic field in the ejecta, and showed that, in 
either case, a bright optical emission (mv ~ 9) is obtained 
at the end of the burst. Meszaros & Rees (1999) extended 
their previous calculations of the RS emission to a radiative 
evolution of the fireball Lorentz factor and pointed out the 
importance of spectral information in constraining the RS 
dynamics and the magnetic field origin from the observed 
t~^ power-law decay of the very early optical light-curve of 
the afterglow 990123 (Akerlof et al. 1999). They also pointed 



out the possibility that optical flashes arise in the same in- 
ternal shocks which generate the burst emission. 

Sari & Piran (1999) have shown that, if the peak fre- 
quency of the RS emission is assumed to be in the optical at 
the time when the optical emission of the afterglow 990123 
peaks (50 seconds after the burst), then the expected soft- 
ening of the RS emission and self-absorption effects can ac- 
count for the radio flare reported by Kulkarni et al. (1999a). 
Kobayashi & Sari (2000) conflrm the RS interpretation of 
this radio flare through numerical calculations of the RS dy- 
namics. 

Chevalier & Li (2000) have presented calculations of the 
RS synchrotron emission until it crosses the GRB ejecta, for 
the case of a wind-like CBM. For their choice of a high mag- 
netic fleld parameter, the RS cooling frequency falls well be- 
low the optical domain, which leads to a RS optical emission 
much dimmer than that observed for the afterglow 990123 
at its peak (40 seconds after the burst). Furthermore, such 
a low cooling frequency implies that the early afterglow op- 
tical emission should cease when the RS has crossed the 
ejecta shell, i.e. at the peak time of the RS emission. Since 
this is in contradiction with the observations of the after- 
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glow 990123, Chevalier & Li (2000) have concluded that a 
wind-like CBM cannot explain the early optical emission of 
the afterglow 990123. 

Constraints on the fireball initial Lorentz factor have 
been obtained by Soderberg & Ramirez-Ruiz (2003) for sev- 
eral afterglows by comparing the observed radio emission at 
~ 1 day with the model RS emission, under the assumption 
that the RS magnetic field and typical electron energy pa- 
rameters (which we shall call microphysical parameters) 
are those determined for the FS from fits to the broadband 
emission of those afterglows (Panaitescu & Kumar 2001). 
Zhang, Kobayashi & Meszaros (2003) have noted that the 
ratios of the RS and FS peak fluxes, and peak and cooling 
frequencies depend only on the fireball initial Lorentz factor 
and the ratio of the magnetic fields, to develop a method of 
constraining these two quantities, and have shown that the 
optical emission of the afterglow 990123 requires a magne- 
tized outfiow. 

In this work we use the general properties (flux, epochs 
during which power-law decays are observed, decay slopes - 
see Table 1) of the radio and optical emissions of the after- 
glows 990123 and 021211 to constrain the ejecta (isotropic- 
equivalent) kinetic energy, CBM density, and the microphys- 
ical parameters for the reverse-forward shock scenario (§4.1), 
for either a homogeneous or wind-like profile of the CBM. 
In contrast with other works, we take into account all con- 
straints arising from the radio and optical measurements of 
the GRB afterglows 990123 and 021211 and we do not as- 
sume certain values for any of the model parameters. We 
also investigate a "wind-bubble" scenario (§5), where all the 
radio and optical emission arises in the FS, with a negligible 
contribution from the RS (which is verified numerically), the 
mitigation of the optical decays observed in the afterglows 
990123 and 021211 being due to the FS crossing the bub- 
ble termination shock, i.e. transiting from a wind-like CBM 
structure to a homogeneous region of shocked wind. For both 
scenarios, we consider either adiabatic or radiative dynam- 
ics, the resulting microphysical parameters being checked for 
consistency with the assumed dynamical regime. 



2 PHYSICAL PARAMETERS AT THE EJECTA 
SHOCK-CROSSING RADIUS 

We begin by calculating the spectral properties (break fre- 
quencies and peak flux) of the RS emission at the radius _/?+ 
where the RS finishes crossing the ejecta shell and the in- 
jection of fresh electrons by the RS ceases. Because most of 
the ejecta kinetic energy has been transferred to the forward 
shock at R+, the radius R+ marks the onset of a steeper 
power-law decrease of the bulk Lorentz factor F with ra- 
dius. After R+, the spectral properties of the RS emission 
can be calculated from the adiabatic evolution of the elec- 
trons and magnetic field. The spectral properties of the FS 
emission can also be calculated from those at R+ or directly 
from the dynamics of the fireball after R+ (i.e. without pass- 
ing through the parameters at _/?+) if the shock dynamics is 
adiabatic. 

Each shock compresses the fluid ahead of it by a factor 
4r' -|-3, where F' is the Lorentz factor of the shocked fluid as 
measured in the frame of the yet unshocked gas, and heats it 
to a energy per particle equal to F' — 1. Therefore, the pres- 



sure equality at the contact discontinuity which separates 
the shocked ejecta and CBM, implies that 



(4F' + 3)(F' - l)nej = (4F + 3)(F - l)n 



(1) 



where F is the Lorentz factor of the shocked fluid in the 
laboratory frame, and n^j, n are the proton number densi- 
ties of the unshocked ejecta and of the CBM, respectively, 
each measured in the corresponding comoving frame. From 
addition of velocities in special relativity. 



r'^FoF(l-,o/3)4(i- + ^) 



(2) 



where Fo is the initial Lorentz factor of the ejecta, /3 de- 
notes velocities and Fq ^ 1, F 3> 1 were used in the ap- 
proximation. Substituting equation (2) in (1), one obtains a 
fourth-degree equation for F, which can be cast in the form 



(Fo - F)^ 
FFo 



(F„ + F)^ 
FFo 



(3) 



Because F < Fo, the first term in the square bracket is at 
least eight times larger than the last term. Ignoring the 1/2 
term, the solution of equation (3) is 



F = 



Fo 



1 + 2Fo 



1/2 



1/2 



(4) 



The limiting cases for the Lorentz factor are 



Fo 



1/2 



/ nej 



1/4 



-^^ > 41 



^ « 4Fi^ 



(5) 



Therefore, in the limit of very dense ejecta, the Lorentz fac- 
tor of the shocked fluid is the same as that of the unshocked 
ejecta, while for more tenuous ejecta, F depends on the ratio 
of the comoving densities. Note that the ratio Uej and ricmb 
changes with the fireball radius R: 



E 



47rmpc2ro(FoA)i?2 



AR' 



(6) 



where E is the fireball ejecta energy (or, if the outflow is col- 
limated, its isotropic equivalent), rrip is the proton's mass, 
A <C .R is the ejecta geometrical thickness measured in the 
laboratory frame (thus FoA is the comoving frame thick- 
ness) , and we restricted our calculations to two simple radial 
structures of the CBM, either homogeneous (s = 0) of parti- 
cle density n or a R~^ stratification (s — 2) corresponding to 
a wind expelled by the GRB progenitor. In the latter case, 
A = 3 X Wr'^ At cm~^, where A* is the mass-loss rate to 
wind speed ratio, normalized to 10~^Mq yr~^/10^kms~^. 
By denoting 



X = 



E 



ATrAmpC^ 



the ratio of the comoving densities is 

Uei X 



AF2i?2- 



(7) 



(8) 



The ejecta-shell thickness. A, is the largest of its initial 
thickness ct, where r is the laboratory frame duration of the 
GRB ejecta release by their source, and the expansion due to 
a spread in the radial outflow velocity of the ejecta particles. 
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This velocity spread can be eitlier a relic of the initial, super- 
Eddington radiation pressure in the fireball, or to an imper- 
fect collimation in the radial direction of the ejecta particles 
at the end of the fireball acceleration. The former leads to a 
comoving-frame expansion of the shell at the sound speed, 
the latter is expected to produce a spread of order l/Fo in 
the ejecta particles direction of motion (Meszaros , Laguna 
& Rees 1993). In either case, the resulting contribution to 
the shell thickness evolution is it/Fg, therefore 




R < 2TlcT 
R > 2rgcr 



(9) 



As shown by Kumar & Panaitescu (2003), the difference 
between the laboratory frame speeds of the unshocked ejecta 
and that of the RS is 



/3o - Prs 



L4 /rgn 

^ 



1/2 



(10) 



for a wide range of the ratio Fon/nej. From here, one can 
calculate the radius R+ at which the RS finishes crossing 
the ejecta shell: 



A{R+) = 



R+ 



Wo - Prs) dR 



(11) 



Once R+ is known, equations (4), (8) and (9) give the 
Lorentz factor F at the shock-crossing radius R+ and all 
the properties of the RS and FS emissions at i?+, which 
can be then extrapolated aX R > R+. We proceed by con- 
sidering separately ejecta shells for which 2rocr < R+ and 
R+ < 2FoCr. In the former case, the ejecta shell undergoes a 
significant spreading while the RS propagates into it, while 
in the latter case the spreading is negligible. The usual ter- 
minology (Sari & Piran 1995) is that of "thin ejecta" for the 
former and "thick ejecta" for the latter. 



2.1 Thin Ejecta Shell - A = R/{2Tl) 

The substitution of equations (8) and (10) in (11) leads to 

1/3 / m \ 1/3 

cm(12) 



(s = 0) : i?+ = ( -!^ ) =0.94x10'^ 



-E53 

fJoi 0,2 



/ r,^ D 0.57X 15 -Ess 

[s — 2) : R+ = — ^yi — = 1.0 X 10 TT- cm 



p2 

t 



A,Tl^2 



(13) 



where the usual notation Qn = Q/10" has been used. The 
defining condition for this case, i?+ > 2Fo cr, becomes 



Fo< 



0.20X 



1/8 



670 



-E53 
nor| 



1/8 



for s — Q 



E53_ 



1/4 



for s 



-<(^)""-»(l 

From equation (8), the density ratio at i?+ is 



^<«*'^]^ 



1.3 Ff, 







3.5 Ff, s = 2 



(14) 



(15) 



(16) 



Comparing with equation (4), this shows that at i?+ the 
Lorentz factor of shocked gas is in an intermediate regime: 



F+ "= 0.60 Fo 



F+ "=' 0.70 Fo 



(17) 



Besides F(J?+), two other quantities are of interest at 
the ejecta shock-crossing radius: the energy of the swept- 
up CBM, Ecbm{R+), and the corresponding observer frame 
time, t+ . The total energy of the shocked CBM is F^ larger 
than its rest-mass energy, thus 



Ecbm(R) = rupc / ^TirAr °F^(r) dr 



(18) 



The arrival time t corresponding to the contact discontinuity 
and the fiuid moving toward the observer is given by 



tcD{R) = (1 + 2) 



dr 
2cF2 



(19) 



where z is the burst redshift. Substituting F from equation 
(4) in equations (18) and (19), one obtains 



(s = Q): Ecbm(R+) = 0.25 S, t+ = 350 (1+^) 



^53 

noFg, 



1/3 
s(20) 



£5 



(s = 2) :£e6m(i?+)=0.34S, i+ = 2.9(l + z)-— ^s.(21) 

/l,i 2 

Thus less than half of the initial ejecta energy has been 
dissipated by the FS by the time when the RS crosses the 
ejecta shell. This means that the shock crossing radius R+ 
is slightly smaller than the usual deceleration radius Rd de- 
fined by E,bm{Rd) = 0.5E. 



2.2 Thick Ejecta Shell - A = cr 

Once again, using equations (8) and (10) in (11), one obtains 
that 



i?4 



{2.0Xcry/^ = 0.24 x 10^^ (^^11^* 
\ no J 



R+ '^^ (0.51Xcr 



)^/=^=0.52xl0i^(:^)'^' 



no 
15 / -E53T0 ' 



(22) 



(23) 



The requirement that 7?+ < 2Fo cr leads to the reversed 
inequalities given in equations (14) and (15). The density 
ratio (equation 8) at R+ is 






(24) 



Because in the thick shell case the ejecta density is lower 
than for a thin shell, the values given in the rhs of equa- 
tion (16) are upper limits for the density ratio. Substituting 
equation (24) in the second regime given in equation (5) 
leads to 



(s = 0) : F+ = 



X 



32(cr)3 



1/8 



= 530 



-E53 
norif 



1/8 



(s = 2) : F+ = 



V4cr/ 



1/4 



= 110 



E. 



A*ro 



1/4 



(25) 



(26) 



Note from equations (25) and (26) that, for a thick 
ejecta shell, the Lorentz factor F of the shocked ejecta is 
independent of that of the unshocked ejecta, Fq. Thus, for 
a sufficiently tenuous ejecta, the contrast between Fo and F 
can be sufficiently large that the RS is relativistic (equation 
2) in the frame of the incoming ejecta. In contrast, in the 
thin ejecta case, equations (2) and (17) lead to F' — 1 = 0.13 
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for s = and T' - 1 = 0.064 for s = 2, i.e. the RS propagat- 
ing in a thin ejecta shell is trans-relativistic. 

The energy dissipated by the FS at i? = R+ and the 
ejecta shell shock-crossing time are 



Ecbm{R+) ^ 0.36 E , t+ = 0.71{l + z)T 



(27) 



for either type of medium. Note that, as for a thin ejecta 
shell, the shock-crossing radius 7?+ is close to the usual de- 
celeration radius. Furthermore, for a thick ejecta shell, the 
observer frame shock-crossing time f+ is fairly close to the 
laboratory frame duration of the ejecta release. Given that 
t+ is roughly the timescale for dissipating the ejecta kinetic 
energy and that, most likely, the shock-accelerated electrons 
cool faster than the dynamical timescale, we expect that the 
duration of an external shock GRB is close to t+. The sim- 
ple temporal structure of the GRBs 990123 and 021211 may 
suggest that they originate in an external shock, neverthe- 
less it is entirely possible that both bursts were produced 
in internal shocks occurring in an outflow with a fluctuating 
ejection Lorentz factor (Rees & Meszaros 1994). In this case, 
the internal shocks take place at the radius Ris ~ V'^i^cr, 
where Fmin is the Lorentz factor of the slower shells; the 
shocked fluid moves at F ~ (FminFmai)^ , where Tmax is 
the Lorentz factor of the faster shells; therefore the observed 
burst duration is t-y — Ris/icF^) ~ {Tmin/Fmax)T. Since a 
high dissipation cfBciency requires a large contrast between 
the Lorentz factors of various ejecta shells (Fmin ^ Fmax), 
the burst duration sets only a lower limit on the duration 
of the ejecta release, r. We shall use this constraint when 
choosing the shock-crossing time t+ for the GRBs 990123 
and 021211 for the case of a thick ejecta shell. 



3 SYNCHROTRON EMISSION 

The synchrotron emission from either shock at any observing 
frequency is determined by the peak flux Fp, the three break 
frequencies - absorption i^a, injection Ui, and cooling Vc - 
and the slope of the afterglow spectrum between its break 
frequencies. 

The peak of the F^ synchrotron spectrum, which is at 
the frequency i/p — mm{i'i, Vc}, is given by 



F^ 



1 



A-Kdl{z 



-N.VPL 



(28) 



where z is the burst redshift, d^ the luminosity distance 
(in a Ho — 70kms~^/Mpc, Q.m = 0.3, f^A = 0.7 universe), 
A^e is the number of radiating electrons, and the factor F 
accounts for the average relativistic boost of the comoving 
frame synchrotron power P^ per electron at the peak fre- 
quency p'p = fp/F. The Doppler boost F appears at only 
the first power in equation (28) because the observer re- 
ceives emission for an area subtending an angle of F^^ radi- 
ans, i.e. from a fraction F~^ of the total number of electrons 
A^e. Since we will be interested in the early time afterglow 
emission, the possible coUimation of the outflow is not an 
issue here. After the shock crossing radius R+, the number 
of electrons energized by each shock is 



N, 



(RS) 



E 



N. 



(FS) 



^" -AR'-' 



3-s 



(29) 



p; = V3(l)p-^B 



P 9 -■ 



(30) 



where cjf>p is the order-unity coefficient calculated by Wijers 
& Galama (1999), e is the electron charge, and B is the 
magnetic field strength. The magnetic field is parameterized 
by the fraction sb of the post-shock energy density stored 
in it. Taking into account that the FS compresses the CBM 
by a factor 4F and heats it is to an energy per proton of 
TmpC? , the magnetic field is 

B = (327r EB AR-' rripc^ F^) ^'^ . (31) 

The break frequencies are calculated from the corre- 
sponding electron Lorentz factors 7a,i,c 

^^^ "" ' BY (32) 



■^p ^ 2 

47r rrieC 



where Xp is another order-unity factor calculated by Wijers 
& Galama (1999) and the last factor in the rhs is for the 
average relativistic boost of the fireball emission by its rela- 
tivistic expansion. The typical electron Lorentz factor after 
shock-acceleration is parameterized as 



7i 






1) 



(33) 



where F' = F for the FS, while F'(i?+) for the RS is that 
given in equation (2) for F = F+ . We assume that each shock 
injects in the downstream fluid electrons with a power-law 
energy distribution 

dN , , _„ 

^(7>70«7 '^ ■ (34) 

The acceleration of new electrons by the RS ceases at R+, 
when all the ejecta have been swept-up, but continues at the 
FS. The cooling electron Lorentz factor is that for which 
the radiative losses timescale is equal to the dynamical 
timescale: 



7c — 67r- 



(35) 



ar (y + l)i?B2 

where ar is the Thomson cross-section for electron scatter- 
ing and Y is the Compton parameter, i.e. the ratio of inverse 
Compton to synchrotron power. The Compton parameter is 
calculated from the electron distribution, as described by 
Panaitescu & Kumar (2000). The random Lorentz factor of 
the electrons radiating at the self-absorption frequency Va is 
given by 



3/10 

'y — / '''P'^P 
'" ^ 1 l/(p+4) 

IP'P 



> 1 



5e Ne 
By^ 47r7?2 



(36) 



The comoving power per electron of equation (28) is 



where 7p = min{7i, 7c} and A^e is given in equation (29). 

Equations (28) - (36) provide the characteristics of the 
RS and FS synchrotron emissions at i? > i?-|-, with the val- 
ues of B, 7i and 7c for the RS calculated at i? = R+ from 
equations (31), (33) and (35), respectively (using equations 
12, 13, 17, 22, 23, 25, 26) and then extrapolated at 7? > _R+ 
as described in §3.2. Once T{R > 7?+) is known, the fireball 
radius can be related with the observer-frame arrival-time 
of the photons emitted at that radius and the evolution of 
the RS and FS spectral characteristics can be calculated. 

The slopes of the piece- wise synchrotron spectrum from 
the broken power-law electron distribution with energy re- 
sulting from shock-acceleration and radiative cooling are de- 
scribed in detail by Sari, Piran & Narayan (1998). Note that. 
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for R > R+ , the RS electron distribution has a sharp cut-off 
at 7c, as the injection of electrons stops when the RS has 
crossed the ejecta shell, which leads to an abrupt switch-off 
of the RS emission at a given frequency when Vc drops below 
that frequency. 



within the first day after the burst. Note that for such high 
radiative losses to occur in an afterglow, FS electrons should 
acquire a substantial fraction of the energy dissipated by 
the shock and should radiate it quicker that the dynamical 
timescale, which requires a sufficiently high magnetic field, 
i.e. a sufficiently large parameter eb and dense CBM. 



3.1 Forward Shock 



If the dynamics of the FS is adiabatic, the Lorentz factor 
Tf of the FS follows immediately from energy conservation 
T%McbmC? = E, where Mcbm is the mass of the swept-up 
CBM. From here one obtains the Blandford-McKee solution 
(Blandford & McKee 1976): 

r^(i?>i?+)=r+f-i^j . (37) 

Equation (19) also gives the photon arrival time for the emis- 
sion arising from a patch on the FS moving at an angle F^^ 
relative to the fireball center - observer axis, from where 
most of the emission arises. Its integration with F from equa- 
tion (37) leads to 



tF{R) 



1 + z R+ 

4-scr! 



_R 

R+ 



+'-! 



(38) 



where the weak deceleration at i? < R+ has been ignored 
in the calculation of the last term in the rhs. For t ^ t+, 
equations (37) and (38) lead to 



rf(f)oct-(3-^'/(«-^»), i?(f) oc t^/(*-^' 



(39) 



Substituting in equations (29), (31), (33), (35) and (36), one 
obtains the following scalings 



z.cOct-(^-'^^'/(«-2^', ^„ OC t-«=/(^«-^») 



(40) 
(41) 



where Va < Vi < Vc was assumed for the last equation. These 
scalings are used to calculate the characteristics of the FS 
synchrotron emission at any f > t+ from those at t = t+. 

If the swept-up CBM radiates half of its internal energy 
faster than the dynamical timescale, the dynamics of the 
afterglow is described by TpM = const, which leads to Vf oc 
jj-(3-3) 'pjjgj-gfoj-e the two extreme regimes of the fireball 
dynamics are 



pk"""; Qj, ^-(3-s)/(8-2s) 



-i(adfa) 
F 



p(™d) j^^-(3-s)/(7-2s) 



(42) 



where the former is the adiabatic case and the latter is a 
highly radiative regime. To estimate the effect of high radia- 
tive losses, we calculate the dependence of the observed flux 
on the fireball Lorentz factor, F^ ex F"^ (with x frequency- 
dependent), and adjust the fluxes obtained in the adiabatic 
case by a factor 

r^(rad) / ^{rad) \ ^ / ^ -(3-s)x/[(8-2s)(7-2s)] 

u 



K 



(adb) 



{t) = 



-^{adb) 



.(43) 



The radiative correction factors at t ~ 0.1 day are close to 
those resulting from the expressions of the afterglow flux 
in the adiabatic case (e.g. equations B1-B9 and C1-C9 in 
Panaitescu & Kumar 2000) if the fireball energy is decreased 
by a factor around 10. This means that a highly radiative 
regime corresponds to a fractional energy loss of about 90% 



3.2 Ejecta Shell 

We calculate the dynamics of the shocked ejecta by as- 
suming adiabatic dynamics and that the ejecta is in equi- 
librium pressure with the energized CBM, whose radial 
profile (Lorentz factor, density, pressure) is described by 
the Blandford-McKee solution. For adiabatic dynamics, the 
pressure in the ejecta shell evolves as 

PR OC {R^'A')-'' (44) 

where A' is the comoving thickness of the ejecta and a is 
the adiabatic index. From the Blandford-McKee solution, 
the pressure equilibrium at the contact discontinuity implies 
that 



PR = Pf{xr) = Pf{x = 1)X_R 



(17-4s)/(12-3s) 



(45) 



where Pf{xr) is the pressure of the shocked CBM at the 
coordinate xr. of the RS, and 



Pf{x = 1) oc Vl-n{R) oc R ^ 



(46) 



is the pressure immediately behind the FS, the last relation 
resulting from the dynamics of the FS (equation 37). Equa- 
tions (44)-(46) imply that the Blandford-McKee coordinate 
for the RS satisfies 



XR OC (i?2-3^/a)3(4-s)/(17-4.) 



(47) 



Equations (37) and (47) and the Blandford-McKee solution 
for the post-FS Lorentz factor 

r(x) « Fpx"'^' (48) 

lead to that the Lorentz factor at the location of the contact 
discontinuity, i.e. the RS Lorentz factor, evolves as 

r_R oc ^-(3-s)/2(-^2a-3^/a^-1.5(4-s)/(17-4s) 



r 



(49) 



The adiabatic index a (equation 49) of the ejecta is ini- 
tially 4/3 if the RS is relativistic, which corresponds to the 
thick ejecta case discussed in §2.2, and lower for a mildly 
relativistic shock (thin ejecta - §2.1), decreasing to 5/3 due 
to the ejecta cooling. Between these limiting cases, the ex- 
ponents of 7? and A' given in equation (49) change by 0.24 
and 0.14, respectively. Such variations do not lead to signif- 
icant changes in the solutions presented in §4.1. For ease of 
further calculations, we will use a = 1.5 in equation (49). 

The only uncertainty left in the evolution of F_r is that 
of the ejecta shell comoving thickness A'. This uncertainty 
also affects the adiabatic cooling of the electrons and the 
evolution of the magnetic field in the ejecta. The evolution 
of the electron Lorentz factors ji and 7c at i? > _/?+ is 

(XV ^ ' OC (R A ) ' 



li, 



(50) 



where V' is the comoving frame ejecta volume and an adia- 
batic index ae = 4/3 for the relativistic electrons has been 
used for the last term. Because the ejecta emission switches 
off when the decreasing cooling frequency Vc falls below 
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the observing frequency u, we will search for afterglow pa- 
rameters for which i/c (tmax) > V, where tmax is the latest 
time when the RS emission was (or thought to have been) 
observed. Therefore the electrons radiating at frequency u 
cool mostly adiabatically and the ejecta radiative cooling 
after R+ can be ignored. For the magnetic field, the flux- 
freezing condition yields 

-2 



Bx ex (i?A' 



Bii oc R 



(51) 



where Bx(B||) is the magnetic field perpendicular (parallel) 
to the radial direction of the fireball motion. 

To assess the effect of the uncertainty of the behaviour 
of A' on the ejecta synchrotron emission, we consider two 
extreme cases: A' = const, as could result from the com- 
pression ejecta against the decelerating contact discontinu- 
ity, and A' = R/Tr, corresponding to a comoving- frame 
expansion of the ejecta shell at a speed comparable to the 
speed of light. For the former case, relating the ejecta radius 
with the observer time through R oc T\t, leads to _R ex f^''* 
for s = and R oc t^''^ for s = 2 (just as for the FS). Substi- 
tuting in equation (51) shows that Bx oc R~^ decays slower 
than B|[. Then, equations (28), (32), (36), and (50) yield 

t n\ u^ ,-0.63 ,-0.96 ,-0.61 /rr,\ 

(s — (j) I'p (X t , l/i^c OC i , Ua OC t (52) 

/ r,\ 771 ,-0.75 ,-1.42 ,-0.72 /I;o^ 

(s = 2) Fp (xt , Ui^c OC i , Ua (xt . (53) 

For an ejecta shell spreading law A' — R/Tr, one obtains 
that R oc i^/^ for s = and R oc t^^'^ for s = 2. Equation 
(51) shows that By oc R^^ decays slower than B±, leading 
to 

/ n\ T^ ,-0.67 ,-1.19 ,-0.41 /r:A\ 

[s = 0) Fp (X t , Ui^c OC t , Ua Oct (54) 

/ r,\ 7^ ,-0.80 ,-1.47 ,-0.47 /r-t-\ 

(s = 2) Fp oat , Vi^c. oat , Ua oct . (55) 

From the above scalings, it can be seen that the temporal 
index of the break frequencies changes by ~ 0.2 (0.05 — 0.25) 
for s = (s = 2), while that of the peak flux by 0.05 for 
either type of medium. The solutions presented in §4.1 vary 
little between the above assumed behaviours of A'. 

Equations (54) and (55) are used to calculate the char- 
acteristics of the RS synchrotron emission at f > t+ from 
those at t+. The effect of radiative losses on the RS emis- 
sion is estimated in a similar way as for the FS (equation 
43), by adjusting the fluxes obtained in the adiabatic case by 
a factor which accounts for the faster deceleration of the FS 
due to the radiative losses. For the highly radiative regime, 
the evolution of the RS Lorentz factor r_R is calculated as in 
the adiabatic case (equations 44 and 49) but using the scal- 
ing of the FS Lorentz factor Tf with radius corresponding 
to the radiative dynamics case (equation 42)^. The absorp- 
tion of the RS radio emission in the FS is also taken into 
account . 



4 REVERSE-FORWARD SHOCK SCENARIO 

The formalism presented in §2 and §3 allows the calculation 
of the RS and FS emission at a given observer time and ob- 
serving frequency. These emissions depend on the dynamics 



of the FS and ejecta shell, i.e. on the fireball kinetic energy 
E and the particle density of the CBM (n for a homogeneous 
medium or A, for a wind surrounding a massive star). The 
ejecta-shell shock-crossing time i+ also depends on the fire- 
ball initial Lorentz factor Fq, for thin ejecta, or on the du- 
ration r of the fireball ejection, for thick ejecta. The initial 
Lorentz factor also determines the number of electrons in 
the ejecta and, therefore, the RS emission. Finally, the RS 
and FS emissions depend on the two microphysical param- 
eters Ei and eb which quantify the typical electron energy 
and the magnetic field. Thus, the RS emission is determined 
by five parameters in the thin ejecta case, or six in the oppo- 
site case, while the FS emission depends on four parameters. 
Note that E and n (or At) determine the emission of both 
shocks. 

In this section we determine in the framework of the 
reverse- forward shock scenario the values of the above pa- 
rameters allowed by the radio and optical emissions of the 
GRB afterglows 990123 and 021211, the only two afterglows 
for which an optical emission has been detected at early 
times, ~ 100 seconds after the burst. 

Table 1 lists the properties of the burst, optical, and 
radio emissions of the two afterglows. For the afterglow 
021211, the optical emission is decaying since the first mea- 
surement, at ti = 130 seconds after the burst (Li et al. 2003). 
For the afterglow 990123, the emission begins to decay at 
~ 45 seconds (Akerlof et al. 1999), after which the burst 
exhibits some variability. This raises the possibility of some 
energy injection in the RS after 45 seconds. For this reason, 
we choose t\ = 73 seconds as the beginning of the after- 
glow decay, as after this epoch the burst exhibits a weaker, 
decaying emission. 

In both cases, the early optical emission fall-offs steeper 
than at later times. For the afterglow 021211, the transition 
between these two regimes has been observed: it occurs at 
t* = 550 - 750 seconds (Li et al. 2003). For GRB 990123, 
the transition is inferred to occur at i* = 400 — 700 seconds 
(Li et al. 2003), i.e. around or after the last early optical 
measurement at f2 = 610 seconds (Akerlof et al. 1999) but 
prior to the next available measurement at ts ~ 4.0 hours 
after the burst (Kulkarni et al. 1999b). 



4.1 Constraints on the FS and RS Emissions 

The slower decaying optical emission lasting for days is natu- 
rally attributed to the FS energizing the CBM. Its t~^ decay 
implies that the FS injection frequency is below the optical 
domain at the first epoch, fs, when the slower decay begins: 



^(ia) < 5 X 10 Hz (constraint 1) 



(56) 



A second constraint is set by the optical flux normalization; 
for this we choose an epoch ti > ta when the model flux, 
Fo , is required to be within a factor 3 of the observed 
flux, Fo " , making allowance for some uncertainty in our 
calculations'^ : 

I Fl;°'"'\t4) < Fi^''^(i4) < 3 F!;°'"'\t4) (constraint 2) .(57) 



^ Because the calculation of F/j makes use of the Blandford- 
McKee solution for adiabatic dynamics, our calculation of F/j 
for radiative dynamics is only a crude approximation 



^ This uncertainty factor determines the width of the region of 
allowed Si — Eb solutions in the lower left corner — upper right 
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Table 1. Properties of optical and radio emissions used to constrain the afterglow parameters 



GRB 


z 
(1) 


(2) 


ti(s) 
(3) 


Fi(mJy) 

(4) 


"12 

(5) 


t2(s) 

(6) 


t3 

(7) 


"34 

(8) 


i4(h) 
(9) 


^4(MJy) 
(10) 


13 
(11) 


t5,6.7(d) 
(12) 


(13) 


990123 
021211 


1.6 
1.0 


70-100 
2-8 


73 
130 


400 
4.1 


1.80 
1.56 


610 
650 


4.2h 
650s 


1.10 
0.94 


8.3 
2.5 


67 
25 


0.68-0.82 
0.55-0.98 


0.2,1.2,4.2 
0.1,0.9,3.9 


130,320,68 
84,44,91 



(1): burst redshift, (2): range of burst duration in various X-ray bands, (3): selected epoch for early optical measurement, (4): i?-band flux 
at ti (Akerlof et al. 1999, Li et al. 2003), (5): ti - t2 temporal index (Li et al. 2003) of the early optical light-curve, Fo oc t-°i2, (6): last 
available measurement or end of the steeply falling-off optical afterglow, (7): beginning of slower decaying optical emission, (8): (3 — t^ 
optical light-curve index (Kulkarni et al. 1999a; Li et al. 2003), Fo oc t~"^i-, (9): selected epoch after (3, (10): i?-band flux at (4 (Fox ct 
al. 2003), (11): slope of optical continuum measured after t-^ (Holland ct al. 2000; Pandey et al. 2003), F^, oc u^l^, (12): selected epochs of 
radio measurements, (13): 20- upper limits on radio fluxes (Kulkarni et al. 1999b; Fox et al. 2003). 



With the exception of the flare seen at t = 1.2 days 
in the radio afterglow of GRB 990123, there are no other 
detections in the radio down to 0.1 mjy or less. We use the 
2(T upper limits on the radio flux at three epochs (^5,6,7) 
spanning the interval 0.1-4 days, to constraint the radio FS 
emission, F, 



(FS), 



F, 



(FS), 



(t5.6,7) < Fy''>{U,6,7) (constraint 3) 



(58) 



That the early afterglow exhibits a steep fall-off requires 
that the RS has crossed the ejecta. The burst duration, t^, 
sets a lower bound on the ejecta crossing-time t+, as dis- 
cussed in §2.2, thus t+ is constrained by 



f^ < t+ < ii (constraint 4) 



(59) 



Furthermore, the steep t^^'^ decay of the RS emission at t\ 
requires that the RS injection frequency is below the optical 
domain at that time: 



(60) 



y'l^^\ti) < 5 X 10^* Hz (constraint 5) . 
Matching the observed flux to within a factor of 3, 

i Fi°''=^(ti) < Fi^^^(fi) < 3Fi°''''(ti) (constraint 6) .(61) 

is another requirement set on the calculated RS emission. 
The detection of RS emission until epoch t2, when the early 
observations end (GRB 990123) or when the transition to 
the FS emission is observed (GRB 021211), implies that the 
RS cooling frequency remains above the optical domain until 
at least epoch t2 



vi"^\t2) > 5 X 10^^ Hz (constraint 7) 



(62) 



otherwise, the RS optical emission would exhibit a sharp 
drop when v^ falls below optical. Finally, the radio upper 
limits are imposed on the RS emission as well: 



F, 



(RS) 



it 



5,6,7 



<F, 



(2<t) 



(is 



6,7) 



(constraint 8) 



(63) 



We search for afterglow parameters (Fo, E; n/A, ■,ei,eB) 
that lead to FS and RS emissions satisfying the constraints 
1-3 and 5-8, respectively. For GRB 990123, constraint 4 re- 
quires that t+ ~ 70 s, as for this burst t^ ~ ti. For GRB 
021211, the same constraint allows that 4s < i+ < 130s. In 
the thin ejecta case, t+ determines the ejecta initial Lorentz 
factor, Fo (equations 20 and 21). For thick ejecta, t+ deter- 
mines the duration of the ejecta release, r, (equation 27). 



corner direction in the flgurcs, and docs not affect significantly 
the conclusions that will be drawn. 



Note that in the case of a thick ejecta shell, the reversed 
inequalities given in equations (14) and (15) provide only 
a lower limit on Fq. Nevertheless, Fo remains a relevant 
parameter because its sets the number of electrons in the 
ejecta. 

The index p of the electron energy distribution (equa- 
tion 34) is determined from the exponent a of the opti- 
cal power-law decay, Fi, oc t~" , for each shock. The avail- 
able measurements of the slope (5 of the optical continuum, 
Fi, oc v~ , at i ~ 1 day, constrain the index p through that 
the intrinsic afterglow spectral slope l3o cannot be larger 
than observed, as intrinsic spectra harder than observed 
(/3o < (3) can be attributed to dust reddening in the host 
galaxy. 



4.2 Results 

The search for afterglow parameters is done by choosing var- 
ious combinations of parameters {E, n) (for s = 0) or [E, A-m) 
(for s = 2) and by identifying the regions in the (e^, eb) pa- 
rameter space which satisfy the above constraints. Various 
values of the initial Lorentz factor Fo satisfying constraint 4 
are tried, to maximize the allowed (e^, eb) parameter range. 
For thin ejecta, the [E, njAt,; Ei,EB) parameter space for the 
RS is significantly smaller than for thick ejecta. For brevity, 
we present here only solutions for the latter case. 

Figures 1-3 show the RS and FS solutions (ei,es) for 
various combinations {E,n) or (i5,^*) for which both RS 
and FS solutions exist. For parameters E, n (or A,) different 
by a factor 10 than those shown, the emission from one of the 
shocks, or from both, fails to satisfy the above constraints. 
For all these three figures, the adiabatic index in the shocked 
ejecta (§3.2) was set to a = 1.5 and the comoving thickness 
of the ejecta was assumed to evolve as A' — R/Fr, where 
F_R is the bulk Lorentz factor of the shocked ejecta at i? > 
R+. Taking a = 4/3, as for a relativistic RS, or assuming a 
constant A' lead to a modest change in the solutions shown 
and leave unaltered the conclusions below. To assess the 
effect of radiative losses, the fiuxes obtained in the adiabatic 
case were decreased as expected from the faster deceleration 
of F in the case of radiative dynamics (equation 43). 

For the reverse-forward shock scenario we reach the fol- 
lowing conclusions about the RS and FS parameters: 

(1) GRB 990123, homogeneous CBM (figure 1): RS and 
FS microphysical parameters differ, with Eg <; 100 Eg 
and El < 0.1 £^ .A RS magnetic field larger than that 
behind the FS points to an ejecta which was initially magne- 
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Figure 1. Reverse-forward shock scenario, GRB 990123, homogeneous circumburst medium (CBM): Reverse and forward shock 
microphysical parameters (eijEs) for the minimal electron energy and magnetic field energy, satisfying the constraints described in §4.1. 
An uncertainty of a factor 3 is assumed in the calculated optical fluxes, however no uncertainty factor is assigned to the analytical radio 
fluxes, as they are compared with 2 a observational upper limits. "RS" and "FSb" denote reverse and forward shock solutions, respectively, 
the FS solutions corresponding to a cooling frequency above (blueward of) the optical domain. The solutions shown are for thick ejecta 
(§2.2), with an observer-frame RS ejecta crossing-time t^ = 70 s, and for highly radiative dynamics (§3.1). There are no FS solutions for 
adiabatic dynamics. Each panel specifies the isotropic-equivalent of the blast- wave kinetic energy, E (in ergs), and CBM particle density, n 
(in cm~'^). For values of {E,n) differing by a factor 10 or more than those shown here, the emission of at least one of the shocks becomes 
incompatible with the observations. For comparison, the isotropic equivalent of the gamma-ray output of this burst is Ej = 3 X 10^** ergs 
(Kulkarni et al. 1999b). 
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Figure 2. Reverse-forward shock scenario, GRB 990123, wind-like CBM: Solutions for radiative dynamics and a CBM shaped by the 
wind blown by a massive star. The parameter A* parameterizes the wind particle density n{r) = M / {Anrripr'^ Vw) = 3 X 10'^^yl*r^, where 
M is the mass-loss rate and v^ is the wind velocity, with A, = 1 corresponding to M/v^ = 10~^ Mq yr~^/10^ kms~^. Higher values of 
the At than shown here, such as those expected for Wolf-Rayet stars (A* ~ 1), do not allow RS solutions, for either adiabatic or radiative 
dynamics. Overlapping FS and RS solutions are shown with larger symbols. There arc no RS solutions for adiabatic dynamics. 



tized. FS solutions correspond to a cooling frequency, i/c , 
above the optical domain. The power-law decay of the FS 
optical light-curve, 034, sets the electron index: p — (4034 -I- 
3)/3 = 2.47±0.07. This implies that the slope of the intrinsic 
optical continuum is jSo = (p— l)/2 = 0.74±0.04, consistent 
with the value reported by Holland et al. (2000), at t = 1 — 3 
days after the burst. FS solutions are obtained only for high 
radiative losses. For the FS parameters shown in figure 1, 
the ratio of the (inverse-Compton) cooling to injection elec- 
tron Lorentz factors is 7c/7i ~ 10 (t/0.1 rf)^ , therefore the 
energy given to electrons which cool radiatively (7 > 7c) is 
a fraction f{t) = (7c/7i)^"'' = 0.34 (f/0.1 d)-^'^ of the total 



electron energy. If the injected electron distribution extends 
to arbitrarily high energies, the fraction of the post-shock en- 
ergy in electrons is e^ = [{p — l)/(p — 2)]ei = 3ei € (0.3, 1). 
Then, at t = 0.1 day, the radiative losses over one dynamical 
timescale are a fraction ^rad(0.1 d) = fe^ ~ e^ £ (0.1, 0.3) of 
the FS energy. If the total electron energy does not exceed 
equipartition (te < 0.5), then ^rad(O.ld) ^ 0.15. Therefore, 
for the FS solutions shown in figure 1, the fireball dynamics 
is between the adiabatic and highly radiative regimes. 

The ejecta kinetic energy E ^ 10^^ ergs (if spherical) 
and the ambient medium density n <; 0.1 cm^"' for which the 
RS emission accommodates the early optical emission of the 
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afterglow 990123 are 10 and 50 times larger than their upper 
limits found by us (Panaitescu & Kumar 2001) from multi- 
wavelength afterglow modelling. The difference is caused by 
the inclusion in the current calculations of the early {t < t,) 
optical emission, which cannot be accommodated by the RS 
if the fireball energy and CBM density were those deter- 
mined from modelling the afterglow emission at i > 4 hours, 
as the constraints 6 and 8 in §4.1 cannot be satisfied simul- 
taneously. Conversely, for the {E, n) values for which the RS 
emission accounts for the early optical emission, the FS so- 
lutions identified in this work satisfy the general constraints 
imposed by observations (§4.1) but provide a x^-wise poorer 
fit to the afterglow data after 4 hours than the best fit ob- 
tained numerically with more accurate calculations of radia- 
tive losses and the integration of received emission over the 
photon equal arrival time surface. 

(2) GRB 990123, r"^ CBM (Hgure 2): there are RS and 
FS solutions with the same microphysical parameters. RS 
solutions exist only for radiative dynamics. Denser wind en- 
vironments are found to allow RS solutions if the kinetic 
energy E is increased, but even for E — 10^^ ergs, which is 
30 times larger than the GRB output, the allowed wind den- 
sity is well below that of a Wolf-Rayet (WR) star {At ~ 1). 
FS solutions correspond to Uc above the optical domain, 
therefore the electron index is p = (4a34-|-l)/3 = 1.83±0.08, 
leading to Po ^ {p — l)/2 = 0.42 ± 0.04, i.e. significantly 
harder than any reported measurement. The RS and FS 
solutions shown in figure 2 with £i ~ 10"^ are consistent 
with high radiative losses if electrons reach equipartition 
(te = 0.5), because p < 2 allows most of the electron energy 
to be at 7 > 7c. However, for the solutions with Si ~ 0.1, 
equipartition energies require a cut-off in the accelerated 
electron distribution at 'ycut < "ic- The synchrotron charac- 
teristic frequency for the 7cut electrons is above the optical 
domain (thus the cut-off does not affect the optical after- 
glow), but 7ctit < 7c implies that all injected electrons are 
cooling adiabatically, i.e. the assumption of high radiative 
losses becomes invalid for £i ^ 0.1. 

(3) GRB 021211, homogeneous CBM : in our previous 
study (Kumar & Panaitescu 2003) we have shown that, in 
the thin ejecta shell case, the RS magnetic field parame- 
ter must be 10^ — 10* times larger than for the FS. For a 
thick ejecta shell, corresponding to a RS shell-crossing time 
longer than the burst duration (2-8 s) but shorter the time 
of the first optical measurement (130 s), there are RS and 
FS solutions with the same microphysical parameters if the 
dynamics is radiative. 

(4) GRB 021211, r-'^ CBM (Bgure 3): RS and FS may 
have the same microphysical parameters, for either adia- 
batic or radiative dynamics. There are also solutions with 
the same parameter Si and a RS magnetic field parameter 
larger than for the FS, indicative of a frozen-in magnetic 
field. FS solutions correspond to Uc above the optical do- 
main, leading to p = 1.59 ±0.08 and f3o = 0.29 ±0.04, which 
is 2a below the hardest slope reported (Pandey et al. 2003, 
at t = 20 hours). Because p < 2, the dynamics may be ra- 
diative for Ei < 0.1 and te = 0.5, while a significantly lower 
te would ensure an adiabatic dynamical regime. Note that 
a high E — 10'''* ergs, more than 100 times larger than the 
burst output, is required by a wind density corresponding 
to At = 0.3, i.e. slightly below that of a WR star. From the 
constraints on the FS parameters. Chevalier, Li & Fransson 



(2004) have also concluded that a weak wind with A^, ~ 0.01 
is required for the afterglow 021211, however, in their cal- 
culations, the FS cooling frequency was placed below the 
optical. If we impose the same constraint here, we do not 
find any solutions for the RS microphysical parameters. 



5 WIND-BUBBLE SCENARIO 

The similarity of the decay indices of the optical light- 
curves of the afterglows 990123 and 021211 before and after 
t* ~ 600 may suggest that a single mechanism produces 
the entire optical afterglow emission. Time-varying micro- 
physical parameters, including the slope of the power-law 
electron energy distribution, could cause a change in the 
optical light-curve decay index, however such an explana- 
tion is in contradiction with the consistency seen in many 
afterglows between the optical spectral slope and light-curve 
decay index at times of order 1 day. 

If the FS is the only mechanism producing the detected 
afterglow emission and if microphysical parameters are con- 
stant, then the non-monotonic behaviour of the optical light- 
curves of the afterglows 990123 and 021211 at t, ~ 600 sec- 
onds must be tied with the fireball dynamics. The fireball dy- 
namics is determined by the ejecta initial energy and CBM 
density. A substantial energy injection can mitigate the af- 
terglow dimming rate, however the energy deposition would 
have to last for the entire duration of the slower power-law 
decay, i.e. until at least a few days, and could lead to a too 
bright radio emission from the RS. 

Besides energy injection, a sudden variation in the ra- 
dial profile of the CBM could also alter the afterglow be- 
haviour. Such a variation is suggested by the association of 
GRBs with the death of massive stars, which drive powerful 
winds, and that the modelling of multiwavelength afterglow 
measurements starting a few hours after the burst leads to 
better fits for a homogeneous medium than a wind. This 
discrepancy may be resolved if the afterglow does not arise 
in a freely expanding wind, but in the environment resulting 
from the interaction of the wind with the circumstellar gas 
(Wijers 2001) or with the winds blown by other stars (Scalo 

6 Wheeler 2001). It is then possible that the GRB ejecta 
run into a CBM whose density profile at smaller radii is the 
r^^ expected for a uniform, free wind, and closer to unifor- 
mity at larger distances. Then, if the FS cooling frequency is 
above the optical domain, the optical afterglow light-curve 
index should decrease by Sa = 0.5 when the wind termina- 
tion shock is reached. The resulting index decrease is slightly 
smaller than observed, but it is possible that deviations from 
uniformity of the environment outside the unperturbed wind 
account for the difference. 

One important issue for this wind-bubble scenario is 
under what conditions the wind termination shock is located 
at the radius R, of the afterglow at the time i* when the 
light-curve transits from a steeper to a slower decay. From 
equations (13) and (38), one obtains 



-R. 



1.4 X 10' 



Es 



U 



A* 600 s 1 ± z 



1/2 



(64) 



_R, is higher by a factor of a few for the wind parameters 
A, ~ 0.1 which we find for this scenario, and by an extra 
factor fZ, 10 for the high ejecta kinetic energy obtained for 
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Figure 3. Reverse-forward shock scenario, GRB 021211, wind-like CBM: Solutions for radiative dynamics and t-^- = 100 s. Adiabatic 
dynamics also allow solutions with the same microphysical parameters behind the shocks. There are no RS solutions for denser winds 
(larger A«), or fireball energies E ^ 3 X 10^^ ergs, for either adiabatic or radiative dynamics. The isotropic-equivalent 7-ray output of this 
burst is B^ = 6 X 10^^ ergs CVreeswijk ct al. 2002). 



the afterglow 990123. Thus, we shall find that R^ ~ 0.3 pc 
for the afterglow 990123 and _R, ^ 0.02 pc for the afterglow 
021211. Because the afterglow radius increases as i? ex f^'* 
a.t R > R, and because the slower decay of the afterglows 
990123 and 021211 is seen from t — t, untill t ~ few days, 
the uniform part of the CBM must extend up to at least 
5i?,. 

Castor, McCray & Weaver (1975) have derived the ma- 
jor physical properties of a bubble resulting from the in- 
teraction of stellar winds with the interstellar gas, taking 
into account the cooling and the diflfusion of the interstellar 
gas into the shocked wind. The radius of the wind termi- 
nation shock Rt can be estimated from the equality of the 
wind ram pressure and that inside the bubble, leading to 
Rt = 4M°-5ii°-_3n(7°'^t5'* pc, where M-5 is the mass-loss 
rate in 10^^ Mq yr^^, ^^,,3 is the wind velocity in km/s, no 
is the interstellar gas density in cm~"^, and and fs is the 
duration of the wind measured in 10^ years. From here, the 
ratio of the contact discontinuity radius Red to that of the 
wind termination shock is Rcd/Rt = 2.3 MZ°-^v'^%no'^t°'^ . 
Thus for a GRB occurring in a dense cloud (n > 10^ cm^^), 
the termination shock radius could be at location required 
by the CBM scenario (equation 64) and the shocked wind 
shell could be sufficiently thick. 

However, WR winds do not interact with the interstellar 
medium but with the wind expelled during the red super- 
giant (RSG) phase, which collides with the main-sequence 
phase wind decelerated by the interaction with the inter- 
stellar medium. The numerical hydrodynamical calculations 
of Ramirez-Ruiz et al. (2001) take into account the wind 
history and show that Rt ~ 0.02 pc for n = 1 cm~^ and 
t = 10® years. Such a termination shock radius is suitable 
for the wind-bubble scenario and the afterglow 021211, how- 
ever the wind-bubble size (0.3 pc) shown by Ramirez-Ruiz et 
al. (2001) is surprisingly small, being 100 times less than that 
expected from the analytical results of Castor et al. (1975). 
Chevalier et al. (2004) have considered the possibility that 
a high interstellar pressure may stall the bubble expansion. 
For the external pressure expected in a intense starburst 
region, their numerical simulations lead to a wind shock ter- 



mination radius Rt — 0.4 pc and a contact discontinuity 
located at Red = 'iRt, which is about right for the wind- 
bubble scenario and the afterglow 990123. 

Alternatively, the uniformity of the CBM at _R > _R* 
required by the wind-bubble scenario might arise from a 
sudden increase in the wind speed, leading to a inner shock 
propagating into the incoming wind. The self-similar solu- 
tions derived by Chevalier & Imamura (1983) for colliding 
winds show that a thick, uniform density shell forms behind 
the inner shock if the termination shock moves at less than 
1% of the unshocked wind speed, which requires the fluctua- 
tion in the wind to consist of a decrease in the mass-loss rate 
by a factor 100 and an increase of the wind speed increases 
by a factor 100 or larger. Such a dramatic change in the 
wind properties exceeds that expected at the transition from 
a luminous blue variable (LBV) wind (M ^ 10~^ Mq yr~^, 
Vyj — 200 km s^^; Garcia- Segura, Mac Low, & Langer 1996) 



or a RSG wind (M = 10"* Mq yr" 



< 100 km s" 



Garcia-Segura, Langer & Mac Low 1996) to a Wolf-Rayet 
wind (M = lO"'^M0yr"\ w„ > 1000 kms"^). The inner 
shock speed Vsh = w™ /lOO = lOkms"^ and the location 
of the termination shock required by the afterglows 990123 
and 021211 imply a WR lifetime R,/vsh = 3 x 10* and 
2 X 10^ years, respectively, i.e. much shorter than predicted 
by evolutionary models for such stars. 

Having found some support for the wind-bubble sce- 
nario in the stalled WR wind-bubble model of Cheva- 
lier et al. (2004), and less so in the interaction between 
the RSG/LR-V and WR winds, we proceed with testing it 
against the radio and optical observations of the afterglow 
990123 and 021211. The wind-bubble scenario must satisfy 
the constraints given in §4.1, the first three pertaining to the 
emission from the fireball interacting with the homogeneous 
portion of the CBM, while the last three referring to the 
FS propagating in the r~^ bubble (instead of the "RS", as 
indicated in those equations). In addition, the constraint 



vi^'^\UM) >5 X 10" Hz 



(65) 



must be imposed, to explain the decrease of the afterglow 
dimming rate at t,, when the uniform medium is encoun- 
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Figure 4. Wind-bubble scenario, GRB 990123: Forward shock microphysical parameters shown separately for the wind-like, inner 
region of the CBM (marked "n oc r~^"), corresponding to the early, fast decaying, optical afterglow (i < 650s), and the homogeneous, 
outer part of the CBM (denoted by "n oc r""), where the t > 4 hours, slower falling-off emission arises. These solutions were calculated 
assuming radiative dynamics, which is consistent with the resulting values of the microphysical parameters for n oc r". Each panel indicates 
the fireball kinetic energy E, the wind parameter A* , and the resulting density n (equation 66) of the uniform outer medium. There are no 
solutions for the the n oc r region for denser winds or lower fireball energies. The lack of overlap between the wind and uniform medium 
solutions indicate that this scenario requires the parameters for magnetic field and minimal electron energy behind the FS to vary when 
the wind-bubble termination shock is encountered. 
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Figure 5. Wind-bubble scenario, GRB 021211: Same as in figure 4 but for the afterglow 021211. The allowed regions for the microphysical 
parameters for the wind and uniform portions of the CBM are shown. In the left panel, the assumption of adiabatic dynamics is justified 
by the resulting total electron energy ee = {p — l)ei/{p — 2) = bsi < 0.1 for the required electron distribution index p = 2.25. In the right 
panel, for the homogeneous CBM region, neither dynamical regime is consistent with the microphysical parameters shown, thus the correct 
n oc r" solutions should lie somewhere between the two extreme regimes, most likely not overlapping with the wind (n oc r~'^) solutions. 



tered. We note that, for s = 2, Vc increases in time, while 
for s = it decreases. Thus, if condition (65) is satisfied, 
then Vc is above the optical domain for any t G (ii,i4). 
We also note that this scenario has only four parameters, 
two {E and >!*) for the FS dynamics and two (e^ and eb) 
for the emission. The density of the uniform region of the 
CBM is determined by the compression by a factor 4 of the 
wind density at the location of the termination shock. From 
equation (64), one obtains that this density is 



m.. 



= 6 X 10 



,Ai 



t* 



£53 V600sl + z 



.(66) 



The density jump across the termination shock should lead 
to a brief brightening of the afterglow (Wijers 2001). Such 



a behaviour may have been missed in the afterglow 990123, 
where there is a gap in the optical observations from 10 min- 
utes to 4 hours. It is not seen in the optical measurements 
of the afterglow 021211 at the "break" time of 10 minutes 
determined by Li et al. (2003), instead an optical emission 
brighter by 0.5 magnitudes than the double power-law fit 
used there is seen at 2 hours after the burst. The lack of a 
brightening at the right time in the afterglow 021211 may be 
problematic for the wind-bubble scenario considered here. 

Figures 4 and 5 display separately the FS microphysi- 
cal parameters which accommodate the observations before 
and after i*, for a few combinations of fireball energy E and 
wind parameter A^. Smaller values of the former parame- 
ter or larger values for the latter do not allow the FS to 
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accommodate the t > t» radio and optical emission of the 
afterglows 990123 and 021211. Note that, because the RS 
emission has at least one free parameter (r and To), requir- 
ing that it does not exceed the measured optical fluxes or the 
radio upper limits, does not constrain the FS parameters. 

As shown in figure 4, the parameters £i and £b satisfy- 
ing the observational constraints cannot be constant across 
_Rt for the afterglow 990123, with the electron energy pa- 
rameter increasing by a factor ^ 100 and the magnetic field 
parameter decreasing by a factor <^ 100 when the fireball 
crosses the wind termination shock. A similar conclusion is 
reached for the afterglow 021211 (figure 5). 



6 CONCLUSIONS 

We have investigated two scenarios that can account for the 
behavior of the early optical emission of GRB afterglows 
990123 and 021211, whose light-curves fall-off as f'^-'^^O-'^ 
at t < i, ~ 600 seconds, while the decay at f > t* fol- 
lows i^io±oi, The first scenario is the widely used reverse- 
forward shock scenario, where the fast decay of the early 
optical emission is attributed to the GRB ejecta energized 
by the reverse shock (RS), and the slower decaying phase 
is associated with circumburst medium (CBM) swept-up by 
the forward shock (FS). 

Figure 1 shows that, for a homogeneous medium, the 
reverse-forward shock scenario can accommodate the ra- 
dio and optical measurements of the afterglow 990123 if 
the ejecta magnetic field is ^ 10 times larger than in the 
shocked CBM, which implies an ejecta frozen-in magnetic 
field (Meszaros & Rees 1997, Zhang et al. 2003), and if the 
RS parameter for the typical electron energy is ^ 10 times 
smaller than for the FS. This differences between the micro- 
physical parameters behind the two shocks are too large to 
be explained away by the inaccuracies in the calculations of 
the afterglow emission presented in §3. 

For a wind- like CBM (figure 2), the reverse-forward 
shock scenario can explain the major properties of the ra- 
dio and optical emissions of the afterglow 990123 with the 
same microphysical parameters behind both shocks (figure 
2), but it requires a wind density corresponding to a mass- 
loss rate to speed ratio less than 10~^ Mq yr~^/10^ kms^^ 
(i.e. A* < 0.1). We obtain a similarly tenuous wind also 
from modelling the broadband data of the afterglow 990123 
at f > 0.2 days, but it should be noted that the best fit with 
a wind-like CBM provides a poorer fit to the data than a 
uniform medium. 

Within the framework of the reverse-forward shock sce- 
nario, the same microphysical parameters are obtained for 
the afterglow 021211, for a thick ejecta shell and either a uni- 
form or a wind-like CBM (the latter case is shown in figure 
3). For a uniform CBM, we obtain n ^ 30 cm""^, larger than 
the n ^ 1 cm^'' inferred by us (Kumar & Panaitescu 2003) 
for a thin ejecta shell. For a wind-like CBM, the afterglow 
021211 requires a wind density corresponding to a mass-loss 
rate to speed ratio below 10~® Mq yr~^/10^ kms~^, a result 
similar to that obtained for the afterglow 990123. 

The second scenario considered in this work (§5), that of 
wind-bubble having a inner r~^ wind- like region surrounded 
by a zone of uniform density, is motivated by that the de- 
crease in the dimming rate of the optical afterglows 990123 



and 021211 seen at t, ~ 650 seconds matches fairly well the 
expectations for such a density profile. The required CBM 
structure finds support in the scenario of WR wind-bubbles 
stalled by a high interstellar pressure discussed by Cheva- 
lier et al. (2004). For this scenario to explain the general 
properties of the radio and optical emissions of the after- 
glows 990123 and 021211, the magnetic field and electron 
energy parameters would have to decrease and increase, re- 
spectively, by a factor of about 100 a,t t — t,, when the wind 
termination shock is reached (figures 4 and 5), a contrived 
feature without a physical foundation. We also find that the 
wind-bubble scenario requires winds which are as tenuous 
as those for the reverse-forward shock scenario. 

Thus the reverse-forward shock scenario provides a 
more natural explanation than the wind-bubble scenario for 
the steep early decay of the optical emission of the after- 
glows 990123 and 021211. Given that the GRB ejecta can 
be initially magnetized and that the RS is less relativistic 
than the FS, the microphysical parameters might differ be- 
hind the two shocks. If their equality is required for a simpler 
scenario, with fewer assumptions, then a wind-like CBM is 
favoured by the reverse-forward shock scenario, though a 
problem still exists: the low wind density inferred in each 
case {A, < 0.1), which is similar to that derived by Cheva- 
lier et al. (2004) for the afterglows 020405 and 021211 and 
by Price et al. (2002) for the afterglow 011211. In the sam- 
ple of 64 Galactic WR stars analyzed by Nugis & Lamers 
(2000) there is only one star with A* < 0.1, the majority 
of the other stars having a mass-loss rate M £ (0.5 — 7) x 
lO~^M0yr~\ a wind velocity v^^ G (1000 - 3000)kms"\ 
and A, G (0.5,3). The dependence of the mass-loss rate 
on stellar mass and metallicity inferred by Nugis & Lamers 
(2000) lead to M ~ IQ-'^ {M /Mq^'^Y^-^ Mqyt-'^ for WN 
stars and M ~ lO'^M /Mq)^-'^Y'^ Z Mq yr-'^ for WCs, 
which may suggest that the tenuous winds required by the 
reverse-forward shock scenario for the afterglows 990123 and 
021211 arise WR stars which are less massive and less metal 
rich than Galactic WRs (Wijers 2001, Chevalier et al. 2004). 
If such stars do not exist, then either the microphysical pa- 
rameters must be different behind the RS crossing the GRB 
ejecta and the FS sweeping-up the CBM or the fast declining 
early optical emission of the afterglows 990123 and 021211 is 
not arising in the RS. One possibility is that the early optical 
afterglow emission is produced in internal shocks occurring 
in an unsteady wind (Meszaros & Rees 1999), a scenario 
which was not investigated in this work. 
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